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A pulse of matter waves may dramatically change its shape when traversing an absorbing barrier 
with time-dependent transparency. Here we show that this effect can be utilized for controlled 
manipulation of spatially-localized quantum states. In particular, in the context of atom-optics 
experiments, we explicitly demonstrate how the proposed approach can be used to generate spatially 
shifted, split, squeezed and cooled atomic wave packets. We expect our work to be useful in devising 
new interference experiments with atoms and molecules and, more generally, to enable new ways of 
coherent control of matter waves. 
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I. INTRODUCTION 

The ability to engineer and manipulate the spatial 
wave function of a quantum particle has far-reaching 
applications in many areas of physics. One example is 
the so-called “beam-splitting”, a coherent division of an 
atomic or molecular wave function into two or more non¬ 
overlapping wave packets (WPs), which is indispensable 
in matter-wave interferometry [TJ . Another example is 
generation of “squeezed states”, i.e., WPs with a reduced 
uncertainty in one dynamical variable (at the expense of 
an increased uncertainty in the conjugate variable) that 
are widely used to enhance the precision of quantum mea¬ 
surements [2]. To date, various mechanisms of reshap¬ 
ing matter WPs have been explored, some of which uti¬ 
lize time-dependent harmonic traps [3], quantum holog¬ 
raphy [3], spatially homogeneous laser light with time- 
dependent amplitude in the presence of a spatially in¬ 
homogeneous magnetic field [5], periodic potentials [6], 
time-dependent electric fields in atom chips [7J, finite- 
length attractive optical lattices with a slowly varying 
envelope [5], or chaotic scattering dynamics [5]. Here 
we report an alternative approach based on the principle 
of time-dependent absorption that has versatile applica¬ 
tions. 



FIG. 1. Illustration of matter pulse carving. 

In this paper, we show that a localized quantum WP 
can be efficiently manipulated - spatially shifted, split, 
squeezed and cooled - by making it pass through a time- 
dependent absorbing barrier, i.e., a narrow region of 


space acting as a particle density sink (see Ref. pTOJ for 
mechanisms and treatment of absorption in quantum sys¬ 
tems). The absorbing barrier can be realized in the con¬ 
text of atom-optics with a sheet of laser light intercepting 
the motion of an atomic cloud (see Fig. [T). The radiation 
frequency of the laser can be chosen such that a passing 
atom becomes undetectable due to ionization or a change 
of its internal state, and thus effectively absorbed. Fur¬ 
thermore, the strength of this absorption process can be 
made time-dependent by varying the intensity of the laser 
in accordance with an externally prescribed function of 
time. 

To further exemplify the process of time-dependent ab¬ 
sorption, we consider an atom initially (at t = 0) pre¬ 
pared in a state ’p(°)(:r)|in), where |in) denotes an inter¬ 
nal state of the atom, and T^^a;) is a spatially local¬ 
ized WP representing its center-of-mass motion. As the 
atom traverses a laser light sheet positioned at x = 0, 
the laser may trigger an atomic transition from |in) to 
another state (or to one of several other states) |abs) of 
the atomic spectrum. The probability of the transition is 
nonzero only inside the light sheet, i.e,. in a close vicin¬ 
ity of the point x = 0, and can be made to depend on 
time by externally modifying the intensity of the laser. 
At some time t > 0, the atom will be found in a state 
'p( t )(:r)|in) + J2 V’ibsO^I 8 ^ 13 )) where ^^(x) and ip^ B {x) 

abs 

denote spatial wave functions corresponding to the inter¬ 
nal states | in) and |abs), respectively. In what follows 
below, we will only be concerned with a projection of the 
full atomic state on | in), and will regard the rest of the 
state, V , ibs( a; )l a ^ )S )’ as a P art that has been removed, 

abs 

or “absorbed”, by the light sheet. It is in this sense that 
we will be interested in a transformation of the center-of- 
mass wave function T, from T* 0 )^) to (x), induced 
by the light sheet playing the role of an absorbing barrier. 

The problem of a non-relativistic quantum particle in¬ 
teracting with an absorbing barrier has a long history 
and is a paradigm of the theory of quantum transients 
m- The limit of an instantaneously opening or clos¬ 
ing barrier, first addressed by Moshinsky, was shown to 
cause well-pronounced oscillations of the probability den- 



2 


sity distribution; this effect is known as “diffraction in 
time” (DIT) and bears close mathematical analogy to 
optical diffraction of light at an aperture with straight 
edges [15]. DIT is generally suppressed if the barrier 
switching is not sufficiently abrupt m or if the parti¬ 
cle momentum distribution exhibits incoherent thermal 
broadening m- In a many particle scenario, the inter¬ 
action between particles is also known to suppress DIT 

mi- 


Here we address the motion of a quantum WP in the 
presence of an absorbing barrier that may, during some 
intervals of time, open or close exponentially fast, but yet 
not fast enough to trigger DIT. We show that this sub¬ 
tle regime, being well within experimental reach, allows 
for controlled reshaping of the WP through the process 
of removing (or “carving out”) parts of the probability 
density profile without a side effect of producing diffrac¬ 
tion ripples. Our analysis takes into account finite tem¬ 
perature effects, which are inevitable in any laboratory 
experiment. 

The paper is organized as follows. In Sec.[TTJ we present 
a general framework for an analytical description of the 
motion of a quantum particle in the presence of a time- 
dependent absorbing barrier. In Sec. m we demonstrate 
the possibility of controlled manipulation of a spatial 
wave function of the particle, and, in particular, show 
how the wa ve function can be d isplac ed (Sec. Ill A I, split 
(Sec. Ill B), and squeezed (Sec. Ill C). We make conclud¬ 
ing remarks in Sec. IV Technical details are deferred to 
two Appendixes. 


II. THEORETICAL FRAMEWORK 

In order to facilitate analytical treatment, we consider 
a quantum particle of mass m initially represented by a 
Gaussian WP 

* ( 0 ) (z)=C,* 0,»oW 

_ ^ 2Re(ao) ^ e -a 0 (:r-xo) 2 4- -i mvo (x-x 0 )/n ^ 

Here, Xq and i>o represent respectively the mean posi¬ 
tion and velocity of the particle, and the parameter a 0 
is related to the spatial extent a of the WP through 
cr” 2 = 2Re(ao). Hereinafter however we assume ao to 
be strictly positive, so that ag = (2a 2 )- 1 . Also, for con¬ 
creteness, we take Xq < 0, a |xo|, and Vq > 0. In 
other words, the initial WP is assumed to be localized 
on the semi-infinite interval x < 0 and moving towards 
the origin. 

After a time t and in the absence of any external forc¬ 
ing, the state ipa},x 0 ,v 0 would evolve into 

/ OO 

dx' Kq } (x - x')^l XO}Vo (x') 

-OO 

= e i ^, Xt , V0 (x), (2) 


where 




2nihT 


exp 



( 3 ) 


is the free-particle propagator, and the functions a T , x T , 
and 4> T are defined as 


ao 


r 1 + 2iHaoT/m 1 
x T = x 0 + v 0 t , 


mVnT 1 

--tan- 


2 H 



( 4 ) 

( 5 ) 

( 6 ) 


We further imagine that in the course of its motion 
the particle encounters an infinitesimally thin absorb¬ 
ing barrier positioned at x = 0. (In a realistic set¬ 
ting, the width of the barrier is assumed to be much 
smaller than the WP size.) The time-dependent trans¬ 
parency of the barrier is characterized by a real-valued 
aperture function yy, ranging between 0 (representing 
zero transparency/complete absorption) and 1 (complete 
transparency/zero absorption). A propagator, transport¬ 
ing the particle probability amplitude from a source point 
x' < 0 to a point x > 0 on the other side of the barrier 
in time t > 0, can be written as muni 


t 

= f drf (^- - Kt T \x)Kf\x '). 

0 

( 7 ) 

The structure of the propagator K^\x,x') stems 
from a superposition of a continuous family of paths, 
parametrized by time r. The propagation of the par¬ 
ticle along each of these paths consists of three consec¬ 
utive stages: First, the particle moves freely from the 
source point x’ < 0 to the barrier at the origin in time 
r; second, the probability amplitude gets modulated by a 
factor proportional to the transparency of the barrier yy; 
third, the particle travels freely to the observation point 
x > 0 through the remaining time (t — r). The additional 

factor 4 — yr) has the meaning of the average clas¬ 

sical velocity at which the particle crosses the barrier, and 
is essential for correctly weighing relative contributions 
of the interfering paths. 

In fact, the propagator K (t \x,x'), given for x > 0 by 
Eq. ([7]), is an exact solution of the following quantum- 
mechanical problem (see Ref. m for full details). In this 
model (x, x'), transporting a wave function from a 
source point x' < 0 to a point x ^ 0, is set to satisfy (i) 
the time-dependent Schrodinger equation 



H 2 d 2 \ 
2m dx 2 ) 


K^ T \x, x') = 0 


( 8 ) 


for 0 < r < t and both x < 0 and x > 0, (ii) the initial 
condition K^ 0+ ^ (x, x') = S(x — x'), (iii) Dirichlet bound¬ 
ary conditions at x = ±00 for negative imaginary times, 
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i.e, A'(~ l l T l)(±oo, x') = 0 , and (iv) two matching condi¬ 
tions relating the propagator and its spatial derivative at 
x < 0 to those at x > 0 , namely 


K^ T \x,x') 

dK^ T \x,x') 

dx 


x=cn 

tc—0“ 


X=CT 


x=0~ 


-(^--Xr)K^\x-x f )\ x=0 , (9) 
dK ( 0 T \x~ x') 


= -(l-Xr)- 


dx 


21=0 


( 10 ) 


for 0 < t < t. Here, Kq denotes the free-particle prop¬ 
agator defined by Eq. ©>■ The matching conditions © 
and ( 10 ) are a time-dependent quantum-mechanical ver¬ 


sion of the absorbing (“black-screen”) boundary condi¬ 
tions proposed by Kottler in the context of stationary 
wave optics |18| ; in their original time-independent for¬ 
mulation, Kottler boundary conditions can be viewed as 
a mathematical justification of Kirchhoff diffraction the¬ 
ory. 

The wave function transmitted through the barrier at 
time t is given by 


OO 

^(z) = o,v 0 ( x ) = J dx ' K(t \ x i^vlx 0 ,v 0 ( x ') 

— OO 

= j dT t (r^ + 77 ) . 

0 

( 11 ) 


where t 0 = |£ol/H) represents the time needed for the cor¬ 
responding classical particle to reach the barrier. Here 
however we choose to focus on a phase-space represen¬ 
tation of the wave function as provided by the Husimi 
distribution OS 


The Husimi distribution quantifies the overlap between 
the time-evolved state X ial,x 0 ,v 0 and a probe Gaussian 
WP centered in phase space around (x, v) and character¬ 
ized by the spatial dispersion a = 1/ y/2ao. Hereinafter, 
we assume 0 < a <C x\ this implies that we only examine 
the wave function fE , ao,x 0 ,uo( :E ) deep inside the transmis¬ 
sion region. Then, using Eq. we obtain 


W 0 ’ liErW ) = dr — 

\t<x o,x,v\^ao,x 0 ,vo! / 9 

0 


a t - T v 


at 


ol t v 0 

a t 0 


X 




OLQ,X,V 


(0) ^ { J 0 \ X0 ,v 0 


(0) 


(13) 


where t = x/v and the asterisk represents complex conju¬ 
gation. The integrand in Eq. (13), unlike that in Eq. o> 
is free of singular points on the closed interval 0 < r < t. 
This fact makes formula (13) an efficient tool for ana¬ 


lytical and numerical inspection of the part of the WP 
transmitted through the absorbing barrier. 

Non-zero temperature, ubiquitous in any laboratory 
experiment, manifests itself as an incoherent broadening 
of the initial velocity distribution of the particle. In order 
to account for this broadening, we describe the particle 
state in terms of a time-dependent density matrix 

( T ) 

Paa,xo : vo,Av 


= — f — e -(«-vo) 2 /Au 2 |^(r) \ /\J/( r ) 

^J_ x Av 1 “0- 


Xq,V I ’ 


(14) 

where Av quantifies the width of the thermal spread of 
the initial velocity. The corresponding finite-temperature 
Husimi distribution reads 

•jyW (x v) = tr ( iw/ 0 ) )(i // 0) A 

rL ao,XQ,vo,Av u ) 1 yr ao,xo,vo,Av' x ao,x,v /ao,x,v\J 
OO 

= 4 / (15) 


Hi alx 0 ,V 0 ( X ,V) = 


\rao,x,v I ctQ ,xq ,vq / 


( 12 ) 


where tr denotes the trace. Substituting Eqs. (12) and 
(13) into Eq. (15), we obtain 


where 


n 


(*) 

CKq 5 3?0 j^0 )Al) V 


(x,v) = 7 [ dr [ dr'XrXr' [^a o ,s,«( 0 )l V’ils.cW 


,X,V K Ot.Q,X,V,X q,Vq,Av ’ 


(16) 


0 0 


oo 

-(v-v 0 ) 2 /Av 2 ( O-t-rV O-tV \ / at-r'V Ut'V \ , ( T ) L (t') /q\ 

ao,x,v,x 0 ,vo,Av ^ J Av ^ Of + a ]M/l 1 ' ~ + ~ 1 ^ao.xovW V<x 0 ,x o,iA U ) 

—OO 


ai a 


|*o IA 


(17) 


In view of the identity ^ v ^ = 77 + 2th J ^ 0 ^ , the integral in Eq. (JTt]) is Gaussian in nature and can be evaluated 
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analytically. (See Appendix [A] for the calculation and 
exact expression.) 

Below we use the Husimi distributions given by 
Eqs. (12) and (13) and Eqs. (16 1 and ( fl7| ) to analyze sev¬ 
eral scenarios of WP engineering corresponding to spe¬ 
cific forms of the aperture function \ T - 


III. WAVE PACKET ENGINEERING 
A. Spatial shifting 

We first investigate physical effects produced by an ex¬ 
ponentially opening (closing) barrier. Thus, we consider 


Xt = Xo e 7r , 


(18) 


where 7 is the rate of change of the barrier transparency 
and xo is its initial value. While an explicit evaluation of 
the corresponding Husimi distribution is not feasible, a 
significant analytical insight can be gained in an asymp¬ 
totic regime, defined by 


1 « N < vtit a_ 

a ~ 2ct 2A ’ 


(19) 


where A = h/(mv 0 ) is the reduced de Broglie wavelength 
of the particle. The first of the three conditions combined 
in Eq. (19), a <C \xq\, has already been used in deriving 
Eq. (11). The second inequality, 2|x 0 | < Vq t, is needed to 
ensure that, at time t , the particle has “passed” the bar¬ 
rier and the dominant part of the WP is well localized in 
the transmission region; under this condition, the Husimi 
distribution of the transmitted WP is accurately repre¬ 
sented by Eqs. ( |T2| ) and (13). Finally, the third inequal¬ 
ity, A < a 2 /(voi), expresses the semiclassical limit, in 
which time variations of the spatial width of the WP can 
be effectively neglected; this limit is commonly known as 
a “frozen Gaussian” regime [2D] , 


Adopting the semiclassical regime specified by Eq. (19) 
and further assuming that 


M < 


2|z 0 K 


( 20 ) 


The average velocity of the particle however remains un¬ 
affected by the barrier. The analysis also reveals that the 
overall transmission probability is given approximately 

by Xt 0 - 

In order to further explore the effect of the WP 
shift, we compute the pure-state and finite-temperature 
Husimi distributions by numerically evaluating the in¬ 
tegrals in Eqs. © and ( p~6|) respectively. To this end, 


we consider a cloud of ultra-cold atoms characterized by 
m = 86.909 u (the mass of a 87 Rb atom), er = 30 /mi, 
v 0 = 3 mm/s, and Av = 0.1 mrn/s. (A cloud of magnet¬ 
ically levitated 8 'Rb atoms with similar parameters has 
been recently used by Jendrzejewski et al. to experimen¬ 
tally demonstrate coherent backscattering of ultra-cold 
atoms in a disordered potential E2-) The cloud is ini¬ 
tially centered at Xq = —0.15 mm and the propagation 
time is set to be t = 100 ms, implying that Xt = 0.15 mm 
= |too I and to = 50 ms = t/2. (In the experiment in 
Ref. [2T], it was possible to let the atomic cloud evolve 
for as long as 150 ms before performing imaging.) Since 
for the chosen set of parameters A cs 244 nm, the sys¬ 
tem is in the semiclassical regime specified by Eq. (19). 


Furthermore, the value of |y| in all computations below 
is taken not to exceed 225 s -1 , which ensures that the 


restriction given by Eq. (20) is fulfilled. 


(a) 


(b) 


(c) 



FIG. 2. 
100 s _1 
curve) 

7 = 100 s 


Aperture function given by (a) Eq. (22 1 with 7 = 
(solid blue curve) and 7 = —100 s 1 (dashed red 


(b) Eq. (231 with 7 = 100 s 
1 — 


and (c) Eq. (24) with 


In the semiclassical regime, the shape of the transmit¬ 
ted WP depends predominantly on the form of the aper¬ 
ture function \ T hr the vicinity of the time to, at which 
the classical particle crosses the barrier, and is largely 
insensitive to the behavior of % T close to the ends of the 
time interval 0 < r < t. So, in order to increase the 
overall transmission probability we consider the aperture 
function (see Fig.[2](a)) 


we use the method of steepest descent to evaluate the 
pure-state Husimi distribution Ha},x 0 , V0 {x,v) analyti¬ 
cally. (See Appendix [B] for details of the calculation.) 

In particular, we show that Ha},x 0 ,v 0 {x,v) is peaked at \ T = min |e 7 *- T_tl ' ) , lj , t\ 
the phase-space point ( x, v ) = (27 + Ax,vq) with 


^ if 7 > 0 
| if 7 < 0 


( 22 ) 


Ax = — 


7(7 

^0 


( 21 ) 


This means that the transmitted WP appears to be spa¬ 
tially shifted compared to the counterpart free-particle 
WP centered at (x t ,vo)- The shift is proportional to the 
rate of change of the barrier transparency and can be pos¬ 
itive (advanced WP) as well as negative (delayed WP). 


when numerically evaluating the integrals in Eqs. (13) 
and (fl6] ) , instead of the one given by Eq. (18). Fig¬ 
ure 131 shows the corresponding Husimi distributions 

Hao,xo,vo {x, V ) and n^l XQ VaAy {x,v) as functions of po¬ 
sition x and velocity v for different values of 7 . Fig¬ 
ure [3ja-c) represent the pure state case, and Fig. §d- 
f) correspond to the case of a mixed, finite-temperature 
state. The spatial shift of the Husimi distribution is well 
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FIG. 4. Husimi distribution a . 0i „ 0i for a mixed 
state, characterized by the same set of parameters as in 
Fig. j3^d-f), in the presence of an absorbing barrier specified 
by Eq. |23| with (a) 7 = 125 s _1 , (b) 7 = 150 s -1 , (c) 
7 = 175 s^, and (d) 7 = 225 s -1 . 


FIG. 3. Husimi distributions for a hl Rb atom in a 
pure (a-c) and a mixed, finite-temperature (d-f) state, i.e., 


H, 


(t) 


ao^O^O (x,v) and v) respectively. System 

parameters are: m = 86.909 u, a = 30 /j, m, xo = —0.15 mm, 
Vo = 3 mm/s, Av = 0.1 mm/s, and t = 100 ms. The aper¬ 
ture function is defined by Eq. (22). 7 = 0 in (a) and (d), 
7 = 100 s -1 in (b) and (e), and 7 = —100 s _1 in (c) and (f). 


pronounced in the figure, and its numerical value is found 
to be in good agreement with the predictions of Eq. (21), 
i.e., Ax = ±30 [im for 7 = ±100 s -1 


It is interesting 
to observe a slight change of the average velocity of the 
particle in the mixed state case (see Fig. [3je,f)). This 
velocity shift stems from the fact that WPs with differ¬ 
ent average velocities, comprising the mixed state, arrive 
at the barrier at different times and, as a result, are sub¬ 
ject to different values of the transparency function. As 
we show later, this effect can be exploited to reduce the 
phase-space uncertainty of (and effectively cool down) an 
atomic cloud. 


Figure [4] shows the response of a finite-temperature WP, 
characterized by the same set of parameters as above, to 
an absorbing barrier specified by Eq. (23). As the rate 7 


increases, the WP stretches and eventually splits in two 
practically non-overlapping parts. A slight difference be¬ 
tween the average velocities of the two parts has the same 
physical origin as in the shifting scenario (see Fig. §e,f)). 

The absorption-based WP splitting mechanism pre¬ 
sented here may be utilized in designing new types of 
matter-wave interferometers. Indeed, the two WPs pro¬ 
duced by the splitting continue propagating along the 
same path in the coordinate space. If this path traverses 
a region with an external potential that is nonuniform 
in space and time, such as a time-dependent disorder, 
then the two WPs will accumulate different phases in 
the course of their motion, and their subsequent recom¬ 
bination will give rise to an interference pattern. The 
interference pattern can subsequently be used to extract 
information about the potential. 


C. Squeezing and cooling 


B. Splitting 

We now consider a different scenario in which the aper¬ 
ture function \ T in the vicinity of to is given by an equally 
weighted sum of an increasing and a decaying exponen¬ 
tial, e 7T and e~ 1T . As before, in order to increase the 
overall transmission probability, we take % T to be unity 
around the ends of the interval 0 < r < t. Thus, we 
choose (see Fig. |2jb)) 

. / cosh[ 7 (r - t 0 )] 1 \ 

** = mm { cosh( 7(0 /2) ’T < 23) 


Finally, we consider a scenario in which the barrier 
first opens exponentially until the time 1 0 and then closes 
exponentially, so that the aperture function reads (see 

Fig.|| c )) 

Xr = e _7|T_to1 , 7 > 0. (24) 

In this case, the Husimi distribution Vo a»(^i 

of a transmitted WP for 7 > 0 appears to be squeezed 
in the x direction and stretched in the v direction com¬ 
pared to the free-particle case. Figure [5] shows the 7 - 
dependence of position (solid red curve) and velocity 
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FIG. 5. Relative spatial (solid red curve), velocity (solid 
blue curve) and phase-space (dashed green curve) uncertain¬ 
ties, and total transmission probability (dashed-dotted purple 
curve) as functions of the rate 7 for the aperture function de¬ 
fined in Eq. ( |24|) . All system parameters are the same as in 
Fig. §d-f) and Fig. [f] 


(solid blue curve) dispersions of the WP, Sx and Sv re¬ 
spectively, computed with respect to the Husimi distri¬ 
bution: 


*(■) = 


J dxdv(-) 2 H — (/ dxdv{-)T-L 


1/2 


(25) 


The initial WP is characterized by the same set of pa¬ 
rameters as above. 

It is interesting to observe that as 7 grows the decrease 
of Sx occurs at a higher rate than the increase of Sv. For 
instance, at 7 = 150 s -1 the spatial dispersion is reduced 
by over 20 % compared to its value in the absence of a bar¬ 
rier, whereas the corresponding relative increase in the 
velocity dispersion is less than 10%. This means that the 
overall phase-space uncertainty 5x{mSv ) decreases with 
growing 7 . The 7 -dependence of the relative phase-space 
uncertainty is shown by a dashed green curve in Fig. [5] 
For the given set of parameters, Sx(mSv) |^_ 0 ~ 3.12 h, 

whereas 5x(mSv) | 7=150s _ 1 — 2.66 H. (We note that the 
Heisenberg uncertainty principle, with all averages com¬ 
puted with respect to a Husimi distribution function, 
states that Sx(mSv) > h jl9].) In other words, the ve¬ 
locity spread of the WP at a finite 7 is closer to the 
Heisenberg limit Svh = H/Sx than that of the corre¬ 
sponding free-particle WP. This in turn means that the 
moving particle gets effectively cooled down by the ab¬ 
sorbing barrier. The cooling occurs through absorption 
of those components of the mixed state that have the 


largest deviations of the velocity from its average value. 
(The dotted-dashed purple curve in Fig. [ 5 ] shows the 
decay of the overall transmission probability defined as 
f+°°dx dvT-L.) In a sense, the effect is similar in 
nature to that of evaporative cooling [22]. 


IV. CONCLUSION 

In summary, we have demonstrated that a moving WP 
of quantum matter can be flexibly manipulated with 
the help of a thin stationary absorbing barrier whose 
transparency changes in time according to an externally 
prescribed protocol. In particular, the WP transmitted 
through the barrier may be spatially shifted, split in two, 
or squeezed and cooled compared to the corresponding 
WP in free space. The reported effects can be observed 
in a laboratory setting using a cloud of ultra-cold atoms 
akin to that produced in experiments in Ref. |21l and a 
laser light sheet of variable intensity. 

In this paper, being mainly interested in a proof-of- 
principle demonstration of absorption-based WP control, 
we have only considered barrier apertures of relatively 
simple, compact functional forms. In real world situ¬ 
ations however aperture function optimization could be 
used to steer the wave function into a desired target state. 
(See Ref. j7] for an example of the optimization approach 
in the context of control of atomic WPs in atom chips.) 
Other important extensions of the present work would be 
to generalize our theory to the case of interacting parti¬ 
cles and to investigate if there are any new effects pro¬ 
duced by a time-dependent absorbing barrier of a finite 
spatial extent. 

We believe that our findings may become of consider¬ 
able value in areas of physics concerned with matter-wave 
interferometry, quantum control, and quantum metrol¬ 
ogy, as well as facilitate better understanding of effects 
of absorption in quantum systems. 
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Appendix A: Evaluation of Wd 

J * OCQ,X,V,XQ,VQ,Z\V 


Here we derive a closed form e xpr ession for the function 

®*Mx 0 ,voAv defined b y Ec f 10- 

Using Eq. ([ 2 ]), we write 
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, ivi (2a 2 .\ X ^ ( , , 2 .mv, . ,mv 2 i 

rj o ] Xo ,v(0) = ( — 1 exp ( -a T {x 0 +vt) - * — (ar 0 + wr) + * 

/ o 2 \ 1/4 

/ 2a: \ \ / . to \ o „ / . m \ 

\ TT + l 2h) TV - 2 [ a T T + l 7^) x ov - ar x o 


V, TTQo / 


exp 


/2«2\ 1/4 

= - exp [ -* 

\7rao/ 


mra T 2 ,mxoa T 


2hac 


-v — i 


han 


v — a T x 0 


In the last line we have used the identity a T r + = where 

Simiiariy, we have 


^ (r ' } (0) 

T OCQ XQ,V \ ) 


2K0 5 

7ra 0 


x exp I i 


vir'a.%, 2 .mx qCX*, 


2Hai 


-v + 1 - 


o 


ha o 


* 2 
v — a T ,x 0 


Therefore 

e-(’-'o) 2/ A" 2 C IO ,„ (0 ) [<?xo,«(°) 


2i?Q;o 


exp (—Au 2 + Bv — C) , 


, 1 .m(rQ T — r'a*,) 

A= — -tt + *- T ’ 


Av 2 


2 han 


2v 0 .mxo (a T — a*,) 

■D — a O ^ j- 

/\v z han 


c — ^2 ( a-r a *') *S > 

_ a T a*, 

-Lb n 


*0 
* \ 2 


a, 


o 


Also, using 


— = V_ .2fiJxo| = v_ _ i 2hxo we write 
*|* 0 |/d “0 rn a 0 m ’ 


a.'j-'u 


ar 


a 


\x 0 \/v 


at — q-fV Oq-'V 


a \x 0 \/v S \ u t 

= Rv 2 + Sv + T 


at — q~v .2hxoa r aq- 

- i - 1 - v 

a+ m ao 


ot*t- T 'V 2hx 0 a*, a*, 

-t -h l -1- V 

a*r m a 0 


r 


where 


S = 


v / a T a t _ T , a t _ r a T 


a o \ a- 


T _ /at-rV _ . 2 hx 0 a T \ l a*t_ T 'V . 2 hx 0 a* T , 


a; 


m 


m 


Substituting the above expressions into Eq. ( fl7| , we ob¬ 
tain 


+oo 

^MxO'VoAv = J dv ( Rv 2 + Sv + T ) ex P (~ Av2 + Bv-C) 

— OO 

1 /2Ra 0 f T BS + R B 2 R\ (B 2 \ 

~~ An V 7r \A^ 2 + 2 A 3 / 2 + 4A 5 / 2 / 6Xf> ~~ )' 


The last expression can be directly adopted for numerical 


evaluation of the thermal-state Husimi distribution. 


















































Finally, we note that, as expected, the last expression 
respects the identity 

lim $ (t ’ T t r ! ) 

<*0 i x ,v ,xo ,vq ,Av 

QL-fc — t'U \ f CX.-fc — 'j-f'V Qt'j-f'Vo 


a+ 


a to J \ ai a to 

X ^i r o ) ,xo^o(°) [<?xo,«o(°)] 
with to = \xo\/vo, thus recovering 

= H «lxo,vo&v) ■ 


Appendix B: Peak of Ha„,x 0 , V0 (x,v) for \ T = Xoe 7T in 
the semiclassical regime 


Here we provide a derivation of Eq. (211. 


In the semiclassical regime, 1 <C ^ <C with 

A = H/(mvo), we define a small parameter 


ht 

e= „ <C 1. 


(e plays the role of an effective Planck’s constant.) Using 


QL-t 


= ao (l + *fe) 1 = ao + 0(e) for 0 < t < t, we write 


/o 2\^/^ / 

1/>& ) ,xo,vo( 0 ) = eX P + Vqt) 


T 1 ) 


1/4 


O(e) 


9 mv n, . .tuvAt 

- — +)+ '^t 


i / ,9 .mvo, . . mvn,T , 

exp ( —a 0 (^o + U)tT - *—^-(a;o +^ot) + » + O(e) 


/ 2a 0 \ 1/4 

(tJ exp 


-aofo^ - ^ao^o^o J T - cto^o “ ^ 


2 . mxo^o 


Similarly, 


u (r,) (0) 

‘ra 0 ,a;o,vo 

'^ (t -A(o)’ 

r an,/ 


2ar 


1/4 


exp 


-a 0 ^o T ' Z “ (^oioxovo - i-^r J r ~ a 0 x 5 + *- ^ 

_ 2/ ,9 __ mv 2 \ , . miro 

—apf (t — t) — 2o.qXv — 1 —— (r — t) — aox + 1 —-— 

\ 2 h ) h 


0 \ / 2 , -mxovo 


1/4 


/ 2a 0 \ 

(tJ exp 

1/4 


(^) ex P 


-9/ / \9 / _ .mv \ . , x ~9 . mx\ 

-ao v (t —t) — I 2ao^ + J ( r ~ t) ~ a o x ~ l ~ 


Then, 


V'io.S.B (°)] V’io.SS,® (°) ^xo.,0 (°) [^L^o (0) 

~ —- exp [ — (7 (t 2 + t' 2 ) — (Vr — iVj) r 
- (Ur + iVi) t’-W ] , 


where 


U = a 0 (7 2 + Uq) , 

Ur = 2a 0 [(x - vt)v + cc 0 n 0 ] 
m (v 2 — Vq) 


V, = 


2 n 


W = 2ao (x — vt) 2 + ; 


r 


The Husimi distribution now reads 

H £lx 0 ,v 0 {x>v) ~ ^ + 4 ^ J dr J dr'xrXr 

0 0 


ib^ r A(0) r d) yT (0) 'ib^ T ^ (0) 

^ olq , x , v \ / V^ao^jt; V w / ‘ra 0 ,tco,^o x ^ f ao^o^o \ / 

t 2 

a 0 (v + Ho) 2 ^_w J rfrXi _ e -Pr 2 -(V B +iV/)r 


27T 


Taking y T = x 0 e 7T , we get 


<>(*,») =, ®>UtA£ e -»' 


27T 


d T e -UT 2 -(V R -'y+iV I )T 


The evaluation of the last integral substantially simplifies 
if we consider the position x to lie sufficiently close to the 
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point xo + vot and the velocity v to be close to i’o. In this 

case, the main contribution to the integral comes from 

the time interval to + St — a/vg < r < to + 6t + a/vo with 

to = ^ and St = |^. (Indeed, since U ~ 2a 0 xg, V R ~ 

4aoXoVo and Vj ~ 0, the exponent — Ut 2 — (V R — 7 )r 

peaks at r max = ~ ^ + ^2 = *o + The 

width of the peak can be estimated as t / -1 / 2 ~ , 1 „ = 

V 2a ov% 

^.) This interval is contained well inside the integration 
range 0 < r < t, provided that |dt| <C <0 or, equivalently, 

ill«A*. 

<7- 

Then, 


mated by Xo, yielding 

H^v^v) ~ Xo e H 

with 

5= (^-7) 2 -F/ 

2 C7 

It is now straightforward (although tedious) to show that 
the exponent S (and so the Husimi distribution) has a 
local maximum at the phase-space point (i m axdmax) = 
(xo + Vot + Ax, no), where 

Ax =_ = 

2a 0 v 0 Vo 

Indeed, one can verify that 




2tt 


+00 


= Xo- 


e -Ur 2 -(y R -'y+iV I )T 


c*o(v + Vq) 2 ( (Vr — 7) 2 - Vj 


2U 


exp 


V 2U 


- IT ) . 


As we are only interested in the form of the Husimi dis¬ 
tribution in the vicinity of the phase-space point (x, v ) = 
(xq + not,no), the exponential prefactor can be approxi- 


dz 

dx 


(x 


max 7 o'max ) 


&Z 

dv 


= 0, 


(5 


max 1 o'max ) 


and 


clet 


d 2 s 

d 2 S 

dx 2 

dxdv 

[I] 

d 2 S 

dvdx 

dv 2 

d 2 Z 


dx 2 

(•£max 7 


_ >0, 


(•Tmax i^max ) 


= -2a 0 <0 . 
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